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Abstract 
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• We clarify the mass dependence of the effective prepotential in N=2 supersymmetric 

SU(A''c) gauge theories with an arbitrary number Nf < 2Nc of flavors. The resulting 
differential equation for the prepotential extends the equations obtained previously for 
SU(2) and for zero masses. It can be viewed as an exact renormalization group equation 
for the prepotential, with the beta function given by a modular form. We derive an explicit 
formula for this modular form when Nf = 0, and verify the equation to 2-instanton order 
in the weak-coupling regime for arbitrary Nf and N^.. 
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I. INTRODUCTION 



New avenues for the investigation of N=2 supersymmetric gauge theories have recently 
opened up with the Seiberg-Witten proposal [1], which gives the effective action in terms 
of a 1-form dA on Riemann surfaces fibering over the moduli space of vacua. Starting with 
the SU(2) theory [1], a form dA is now available for many other gauge groups [2], with 
matter in the fundamental [3] [4] or in the adjoint representation [5]. This has led to a 
wealth of information about the prepotential, including its expansion up to 2-instanton 
order for asymptotically free theories with classical gauge groups [6] . 

These developments suggest a rich structure for the prepotential JF, which may help 
understand its strong coupling behavior, and clarify its relation with the point particle 
limit of string theories, when gravity is turned off [7] . Of particular interest in this context 
are the non-perturbative differential equations derived by Matone in [8] for SU(2), and 
later extended by Eguchi and Yang in [9] to SU(A^c) theories with only massless matter. It 
was however unclear how these equations would be affected if the hypermultiplets acquire 
non- vanishing masses. 

In the present paper, we address this issue by providing a systematic and general 
framework for incorporating arbitrary masses rrij. In effect, the masses mj are treated 
on an equal footing as the vev's of the scalar field in the chiral multiplet, since they 
are both given by periods of dA around non-trivial cycles. For the masses, the cycles are 
small loops around the poles of dA, while for Ofe, they are non-trivial ^-homology cycles. 
This suggests that the derivatives of with respect to the masses should be given by the 
periods of dA around "dual cycles" , just as the derivatives of with respect to are given 
by the periods of dA around S-cycles. We provide an explicit closed formula for such a 
prepotential, motivated by the r-function of the Whitham hierarchy obtained in [10]. (In 
this connection, we should point out that intriguing similarities between supersymmetric 
gauge theories and Whitham hierarchies had been noted by many authors [11], and had 
been the basis of the considerations in [9], as well as in [4], the starting point of our 
arguments). Written in terms of the derivatives of JF, this closed formula becomes the 
non-perturbative equation for that we seek. It can be verified explicitly to 2-instanton 
order, using the results of [6]. 

Specifically, the differential equation for is of the form 
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[Resp_{zdX)Resp_{z-^dX) + Resp^{zdX)Resp^{z-^dX)] 



(1.1) 



27ri 



with V the operator 




(1.2) 
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The right hand side in (1.1) has been interpreted in [8] [9] in terms of the trace of the 
classical vacuum expectation value X^^i a|, although there are ambiguities with this 
interpretation when Nf > Nc- Mathematically, it can be expressed in terms of ■i?-functions 
for arbitrary Nc when Nf = (c.f. Section III (c) below). There is little doubt that this 
should be the case in general. Now we have by dimensional analysis 

{V + A^)J^ = (1.3) 

if A is the renormalization scale of the theory. Thus the proper interpretation for the 
equation (1.1) is as a renormalization group equation, with the beta function given by a 
modular form! 

Finally, we observe that the effective Lagrangian in the low momentum expansion 
determines the effective prepotential only up to afc-independent terms. However, masses 
can arise as vacuum expectation values of non-dynamical fields, and we would expect 
the natural dependence on masses imposed here to be useful in future developments, for 
example in eventual generalizations to string theories. 



II. A CLOSED FORM FOR THE PREPOTENTIAL 

(a) The geometric set-up for N=2 supersymmetric gauge theories 

We recall the basic set-up for the effective prepotential J- of N=2 supersymmetric 
SU(A^c) gauge theories. 

The moduli space of vacua is an Nc — 1 dimensional variety, which can be parametrized 
classically by the eigenvalues a^, Xl^i Ofc = of the scalar field (p in the adjoint represen- 
tation occurring in the N=2 chiral multiplet. (The flatness of the potential is equivalent 
to [(^, (/>^] = 0). Quantum mechanically, the order parameters cik get renormalized to pa- 
rameters afc. The prepotential JF determines completely the Wilson effective Lagrangian of 
the quantum theory to leading order in the low momentum expansion. Following Seiberg- 
Witten [1], we require that the renormalized order parameters a^, their duals an^k, and 
the prepotential JF be given by 

dA, dD k = S 

27rz /^^ ' 2771 Jb^ 

0'D,k (2.1) 

where dA is a suitably chosen meromorphic 1-form on a fibration of Riemann surfaces V 
above the moduli space of vacua, and Aj^ Bj is a canonical basis of homology cycles on F. 



_ 

dak 
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In the formalism of [4], the form dA is characterized by two meromorphic Abehan 
differentials dQ and dE on F, with dA = QdE. For SU(A^c) gauge theories with Nf 
hypermultiplets in the fundamental representation, Nf < 2Nc, the defining properties of 
dE and dQ arc 

• dE has only simple poles, at points P+, P_, Pi, where its residues are respectively 
— A^c, Nc — Nf, and 1 (1 < i < Nf). Its periods around homology cycles are integer 
multiples of 27rz; 

• Q is a well-defined meromorphic function, which has simple poles at and P_, 
and takes the values Q{Pi) — —mi at Pi, where mj are the bare masses of the Nf 
hypermultiplets ; 

• The form dA is normalized so that 

Resp+(zdA) = -Nc2-^/^% Resp_(zdA) = {N^ - Nf)^^— )i/(Afc-iv,) 

Resp+(dA) = (2.2) 

where A is the dynamically generated scale of the theory, and z = E~^^^'' or z = 
^i/{Nc-Nf) jg ^YiQ holomorphic coordinate system provided by the Abelian integral E, 
depending on whether we are near P^ or near P_. 

It was shown in [4] that these conditions imply that F is hyperelliptic, and admits an 
equation of the form 

v''={\[iQ- ~^k)f - A'^=-^^ + ^ ^(^)' - ^(^) (2.3) 

fe=l 3=1 

Here cik are parameters which coincide with when N^. < Nf, but may otherwise receive 
corrections. It is convenient to set 

The function Q in dA = QdE is now the coordinate Q in the complex plane, lifted to the two 
sheets y — ■^\/ A? — B of (2.3), while the Abelian integral E is given by i? = log {y + A(Q)) . 
The points P± correspond to Q = oo, with the choice of signs y — ±\/ — B. 

(b) The prepotential in closed form 

We shall now exhibit a solution J-' for the equations (2.1) in closed form. Formally, it 
is given by 

- Nf p. 

Uk (p dA — nij / dA 



1 

2T = 

2m 



fc=l "^k j = l 

+ Res (^dA)Res p^ {z~^dX) + Res p_ (2;dA)Res p_ {z~^d\) 



(2.4) 
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However, the above expression involves divergent integrals which must be regularized. For 
this, we need to make a number of choices. First, we fix a canonical homology basis Ai, 
Bi, along which the Riemann surface can be cut out to obtain a domain with boundary 
Yi^J^i^ A~^B~^AiBi. Next, we fix simple paths C_, Cj from to P_, Pj respectively 
(1 < j < Nf), which have only P+ as common point. As usual the cuts are viewed as 
having two edges. With these choices, we can define a single-valued branch of the Abelian 
integral E in Tcut = T \ (C_ U Ci U ■ ■ ■ U Cat^) as follows. Near the function 
provides a biholomorphism of a neighborhood of to a small disk in the complex plane. 
Choose the branch of logQ~^ with a cut along Q~^{C-), and define an integral E of dE 
in a neighborhood of P+ in Tcut by requiring that 

E = NclogQ + log2 + 0{Q-'^) (2.5) 

The Abelian integral E can then uniquely defined on Tcut by integrating along paths. 
It determines in turn a coordinate system z near each of the poles P+, P_, and Pj, 
l<j<Nf,e.g., 

z = e~^^ near P+ (2.6) 

It is easily seen that z is holomorphic around P_|_, and that z = 2^Q~^ + 0((5~^)- The 
next few terms of the expansion of z in terms of Q"^ are actually quite important, but we 
shall evaluate them later. Similarly, we set z = e^^"^/ near P_, and z = near Pj, 
l<j<Nf. 

The same choices above allow us to define at the same time a single-valued branch of 
the Abelian integral A in Tcut- Specifically, A is defined near P_|_ by the normalization 

A(^) = -Resp^ (zdX) - + 0{z) (2.7) 

z 

with z the above holomorphic coordinate (2.6). As before, A is then extended to the whole 
of Tcut by analytic continuation. Evidently, near P_ , A can be expressed as 

A(^) = -Resp_(^dA)- + A(P_) + 0{z) (2.8) 

z 

in the corresponding coordinate z near P_, for a suitable constant A(P_). Similarly, near 
Pj, A can be expressed as 

A(^) = -mj\ogz + X{Pj) + 0{z) (2.9) 

for suitable constants Pj. The expression (2.4) for the prepotential can now be given a 
precise meaning by regularizing as follows the divergent integrals appearing there 

dA = \{Pj) - A(P_) (2.10) 
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This method of regularization has the advantage of commuting with differentiation under 
the integral sign with respect to connections which keep the values of z constant. 



(c) The derivatives of the prepotential 

The main properties of T are the following 





1 


dak 


~ 271^ 


dT 


1 


drrij 


~ 27H 



dX 



1 



(3) 



(2.11) 
(2.12) 



where dO^ are Abelian differentials of the third kind with simple poles and residues +1 
and -1 at P_ and Pi respectively, normalized to have vanishing Aj-periods. We observe 
that the Wilson effective action of the gauge theory is insensitive to modifications of T by 
Ofc-independent terms. The equation (2.12) can be viewed as an additional criterion for 
selecting T ^ motivated by the fact that the mass parameter — of dA can be viewed as a 
contour integral of dA around a cycle surrounding the pole Pj. In analogy with (2.4), the 
derivatives with respect to of a natural choice for T should then reproduce the integral 
of dA around a dual cycle. This is the origin of the first term on the right hand side of 
(2.12), if we view the path from P_ to Pj as such a dual "cycle". The second term on the 
right hand side of (2.12) is a harmless correction due to regularization. The expression 
between parentheses is actually always a multiple of 7ri, although we do not need this fact. 

We now establish (2.11) and (2.12). We need to consider the derivatives of dA with 
respect to both and vrij. We use the connection = V of [4], which differentiates 
along subvarieties where the value of the Abelian integral E (equivalently the coordinate 
z) is kept constant. Then simply by inspecting the derivatives of the singular parts of dA 
in a Laurent expansion in the 2;-coordinate near each pole, we find that 

Va.dA = 27rzda;fe, V^.dA = ddf\ (2.13) 

where dcufc is the basis of Abelian differentials of the first kind dual to the Aj^-cycles. 
Next, we recall from (2.2) that the residues Resp^(2;dA) and Resp_(2;dA) are constant. 
Consequently, 

da;fc - ^ m j / da;fc 

>+(2;~Ma;fc) + Resp_(2;dA)Resp_(2;~Ma;fe) (2.14) 



dak 



1 
27ri 



dA + ^ ai 



i=l 



+ Res p. (2;dA)Resj 
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However, we also have the following Riemann bilinear relations, valid even in presence of 
regularizations 



duji 



1 

1 
27H 



J dojk 



Resp_,_(2; ■'^da;^) 



-Resp^{z-^dnf^) 



(2.15) 



Here dQ)^' are Abelian differentials of the second kind, with a double pole at P±, vanishing 
^-cycles, and normalization 

(2.16) 



dn^l^ = z-^dz + o(z) 



The relations (2.15) follow from the usual Riemann bilinear arguments, by considering 
respectively the (vanishing) integrals on the cut surface Tcut of the 2-forms d(a;jda;fc), 
d{nf^dLOk), d{nf^diOk), d{nf^dn±). Applying (2.15) to (2.14), we obtain 



2-— =— i dA + V tti i doji + — S2mj (p dnf^ 



+ — Resp, (^dA) / doL^^ + — Resp_(^dA) / dVP^ 
27rz Jb^ 27rz Jb^ 



(2.17) 



However, the expression 



Nr. 



d\ = 2Tii ^ ttidui + Resp^ {zdXjdVt^l^ + Resp_ {zdXjd^i:^ + rUjdVL 



)(2) ^ V^_.^q(3) 



i=l 



3 = 1 



(2.18) 



is just the expansion of dA in terms of Abelian differentials of first, second, and third kind! 
The equation (2.11) follows. The equation (2.12) can be established in the same way. First 
we write 



dJ^ 1 



dmi 2TTi 



„ „p, Nf 

^aij> dOf ^ - / dX-J2^j dnf^ 



+ Resp+(2dA)Resp^(2-Mnf ^) + Resp_ (^dA)Resp_ (^"M^^'^O 



(2.19) 
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Substituting in the bilinear relations gives 



dT 1 



dmi lui 



— 27ri 2_] o-i / djJi — dX — rrij / dfl. 



-Resp^{zdX) I dn^^^ -Resp_{zdX) d^L^^ 



Pi 



(2.20) 



Again, the Abelian differentials recombine to produce dA, and the relation (2.12) follows. 



III. THE RENORMALIZATION GROUP EQUATION 

(a) The renormalization group equation in terms of residues 

Combining the equations (2.4), (2.11), and (2.12) gives a first version of the renor- 
malization group equation for J^, valid in presence of arbitrary masses mj 



EOJ^ \ - 



k=l 



1T = 



27rz 



Resp^(2dA)Resp^(2 ^dA) 
+ Res p_ (2;dA)Res p_ {z~^dX) 



(3.1) 



(b) The renormalization group equation in terms of invariant polynomials 

We can evaluate the right hand side of (3.1) explicitly, in terms of the masses m^, and 
the moduli parameters and A of the spectral curve (2.3). For this, we need the first 
three leading coefficients in the expansion of Q in terms of z at and P_. Now recall 
that at P+, (5 ^ oo, y = \l — B, and 

^=(y + ^)-V^= (3.2) 

For Nf < 2Nc, we may expand — B in powers of B /A^ and write, up to 0{Q^''~^) 

(3.3) 



, ^r, IB 1 S2 

y + A = 2[A--^--^] 



We consider first the terms in (3.3) of order up to 0{Q^--'^). Then for Nf < 2Nc - 2, 
only the top two terms in A contribute, while for Nf = 2Nc — 1, we must also incorporate 
the term K^x^f-^- = k^x^-'^ from B/A. Thus 

A + y = 2g^= [1 - (si + 5n„2n^-i^)Q-^] + o(g^=-2) 
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where we have introduced the notation 



fci<---<fci 



This leads to the first two coefficients of z in terms of Q, or equivalently, the first two 
coefficients of Q in terms of z 

Comparing with (2.2), we see that this confirms the value of Resp^(2;dA) required there, 
while the condition that Res (dA) = is equivalent to 



A2 



Nf,2N^-l- 







(3.4) 



Similarly, in the expansion of A + y to order 0(Q^'=~^), we must consider separately the 
cases Nf < 2Nc — 2, Nf — 2Nc — 2, and Nf — 2Nc — 1, depending on whether the terms 
B/A and B'^/A^ contribute to this order. Taking into account (3.4), we find 



Q = 2-1/^=^-1(1 



(3.5) 



with defined to be 



S2 = §2 — Sn^^2N2-2-^ SNf,2N^-l-^ti 



(3.6) 



Near P_ , we have instead 



. B ,1/2 IB IB^ 1 B^ 

A + y = A-Ail-^y/^ = -- + -^ + -- 

Again, considering separately the cases Nf < 2Nc — 2, Nf = 2Nc — 2, Nf — 2A^c — 1, we can 
derive the leading three terms of the expansion of 2; = E~^/^^f~^'=^ = {A + y)~i/(%~-'^c) 
in terms of Q. Written in terms of an expansion of Q in terms of z, the result is 



A^ 



.l/{Nf-N,)-l 



1 - 



Nf-N,' 2 



(A!)V(a^/-JVc), 



+ 



1 



{Nf - N,r ' 2 
with S2 given by 



i^f/^^f-^^\S2iNf - AT,) + 1(-1 + Nf- N,)tl)z' 



(3.7) 



A^ A^ 

S2 — t2 — diV/,2iV2-2— j 0Nf,2Nc-l—rtl 



(3.8) 
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Since dA = 



-NcQ^ near P+ and dA = 



-{Nf -Nc)Q^ near P_, we obtain 



Resp+(^-MA) =2^/^=52+, 

Resp_(.-MA) = - (^)^/(^^-^=)(52- + ^(1 



1 




(3.9) 



Subsituting in the values of Resp^(2dA) and Resp_(2;dA) given in (2.2), and rewriting the 
result in terms of S2 and the operator T> of (1.2), we can rewrite the renormalization group 
equation (3.1) as 



Before proceeding further, we would like to note a few features of the renormalization 
group equation and of our choice of prepotential. 

(1) The RG equations (3.1) and (3.10) are actually invariant under a change of cuts. In- 
deed, a change of cuts would shift the values of the regularized integrals (1.4) by a 
linear expression, and hence by a quadratic expression in the masses rrij , indepen- 
dent of the afc. In view of Euler's relation, such terms cancel in the left hand side of 
(3.1) and (3.10). Thus the right hand side of the RG only transforms under a change 
of homology basis, and is a modular form; 

(2) From the point of view of gauge theories alone, we can in practice ignore on the right 
hand side of (3.1) and (3.10) terms which do not depend on the afc. Such terms can 
always be cancelled by a suitable afc-independent correction to JF. These corrections 
do not affect the Wilson effective action since it depends only on the derivatives of 
with respect to a^; 

(3) Some caution may be necessary in interpreting §2, in terms of the classical order 
parameters ak- In particular, when Nf > Nc, there are several natural ways of 
parametrizing the curve (2.3), which the Ofc get shifted in different ways to Ofc 7^ Ofc 
[3] [4] . As noted in [6] , the prepotential is independent of such redefinitions of the 
ak- However, this would of course not be the case for S2 = 5^^^^ Ofca^, which argues 
for a distinct interpretation for §2 = '^j^k ^kO-j- 

(c) The renormalization group equation in terms of i?-functions 

As noted above, the right hand side of the RG equation (3.1) is in general a modular 
form. For Nf = (and arbitrary Nc), we can exploit the symmetry between the branch 
points given by = (^ — A)(^ + A) = Y[k=i(Q ~ ^k)(Q ~ ^k) known formulae for 



27ri'DJ^ = -{Nf - 2Nc){s2 - SNf,2N^-2^ - (^7Vj,2JVc-l 

+ (Nf-Nc)t2-l{Nf-Nc-l)tl 




(3.10) 
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their cross-ratios to write it explicitly in terms of i^-f unctions. More precisely, we observe 
that 



E4 = B4f = E(4)' (3-11) 



fc=l fc=l k=l 



Let the canonical homology basis be given by cycles surrounding the cut from to 
x'l^ , 1 < k < Nc — 1 on one sheet, and by cycles going from x]^^ to on one sheet, 
and coming back from x^ to on the opposite sheet. Then for the dual basis of Abelian 
differentials da; = (da;fc)fe=i,...,jv^_i, we introduce the basis vectors e^^^ and t^''^ of the 
Jacobian lattice by 

JAk JBk 

We have then the following relations between points in the Jacobian lattice 

du = Je(^), da; = ^(r^^+i) + r^'^)) (3.12) 

Let 0(Q) denote the Abel map 

rQ rQ 
(p{P)^{ dioi,---, diON^-i) 

JQo JQo 

If we choose Qq so that = \t^^\ it follows from (3.12) that 







+ • 






<t>{4) 




+ • 


•• + e('^)) + -r^'^), 1 <A;<iVc-l 








+ • 








=0 






(3.13) 



If we introduce the functions Fi{Q) by 



an inspection of the zeroes shows that we have the following relation between and 
cross-ratios 

Ff(Q') _ Q'-xi Q-x^^ 



X 



(3.15) 



Nr. 
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For the Riemann surface (2.2), we also have for all Q 



HiQ - x+) = A{Q) - A = HiQ -x^)- 2A 



(3.16) 



1=1 



1=1 



Setting Q = x'l^ gives the relation 



Nr. 



n(4-^r) = 2A 



(3.17) 



1=1 



Combining with products of expressions of the form (3.15) evaluated at branch points, we 
can actually identify the branch points 

X+-X+ =A{Gk-Gi) 



A 



Nr. 



1=1 



(3.18) 



where is defined to be 



(3.19) 



Since FI^{xj^ ) is independent of A;, this expression may be simplified. 



Nr Nr 



1=1 k'=i 



(3.20) 



The evaluation of the functions Fl' on the branch points is particularly simple, and we 
have 

i^icPix-) + (P{x-) + (P{x+)\t)^ 



Fnx;,) = 



(3.21) 



^?((/>(x^J+(/)(x^) + (/)(4)|r)2 

where the values of (j){x^) can be read off from (3.13). This leads to the following expression 
for the right hand side of (3.10) : 



Nr 



k=l 



Nr 



A^ 



k=l 



(3.22) 



k=l 



which is a modular form. 
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IV. THE WEAK-COUPLING LIMIT 



It is instructive to verify the renormalization group equation (3.10) in the weak- 
couphng hmit analyzed in [6] to 2-instanton order. 

We recaU the expression obtained in [6] for the prepotential T to two-instanton order 
in the regime of A — > 0. Let the functions Six) and S^ix) be defined by 

^ ng.(x + m,) ^ 

Then the prepotential T is given by 

with the terms .F^^^ T^'^^ corresponding respectively to the one-loop perturbative 

contribution, the 1-instanton contribution, and the 2-instanton contribution 

2.i^(«) = - i Et"- - <")'i<'8 ^^^^^ + i E(«' + ir'' 



fc=i 

2.z^(^) ^lA4(^ ^^K)^KaO ^ \ Y^s,{a,)dlM-^)) (4.2) 
16 (afe-a,)' 4^ 

Here we have ignored quadratic terms in Ofc, since they are automatically annihilated 
by the operator V. We also note that the arguments of [6] only determine T up to a^- 
independent terms, and thus we shall drop all such terms in the subsequent considerations. 
The formulae (4.2) imply 



E + E-.^ - 2^ = - 2iVJ(^ Eal + ^« + 2^<^.) ,4.3) 



where all A^ terms have been ignored. 

On the other hand, up to a^-independent terms, the renormalization group equation 
(3.10) reads 

E«5^ + E'».a;;^-2^=«(^/-2^JE«? (") 

k=\ j=l ■' k=l 
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where we have rewritten §2 as 

No TO 

S2 



1 ^= A2 

- 2 I] ^fc + J^^Nf,2N,-l (4.5) 
A:=l 



To compare (4.3) with (4.4) we need first to evaluate ^^li a'l in terms of the renormahzed 
order parameters ak- Using the formula (3.11) of [6], this can be done routinely 

ak = ak + ^dkSki&k) + ^dlSkiak) + 0(1^) (4.6) 

where we have set dk = d/dak, and defined functions S{x), Sk{x) in analogy with (4.1), 
but with ak replaced by dk- Inverting dk in terms of ak, and rewriting the result in terms 
of the derivatives dk = d/dak with respect to the renormahzed parameters ak, we find 



1=1 



ak = ak- —dkSkiak) - -^dlSk{akf + — ^ az5Kaz)afea^5fc(aO + 0(A6) (4.7) 



and hence 



^2 No ^4 iVc 

= ^al - — akdkSk{ak) " XI <^kdlSk{akf 

k=l k=l k=l k=l 

+ Y E akdiSi{ai)dkdiSk{ak) + — Y,i9kSk{ak)f + 0(A6) (4.8) 

k,l=l fc=l 

Next, we need a number of identities which can be established by contour integrals, in 
analogy with the identities in Appendix B of [6] 

ttkOkSkittk) = — X^ Sk{ak) + {afc— independent terms} 

k=l k=l 

Y,(^kdlSk{akf = -^Y.^lSk{akf 

k=l k=l 

J2akdiSiiai)dkdiSk{ak) = - 2 1^(9,5, (a,))^ + 3^ ^^M^ 

-J2Sk{ak)dlSk{ak) (4.9) 

k^l 

Using (4.9) we can indeed recast ^k=i ^k ^ 
E - E 4 + E + ^ ( E ^t'-a^^ + I Sk{ak)dlSk{ak)) (4.10) 

k=l k=l k=l k^l ^ ^ k=l 

The equality of the two right hand sides in (4.3) and (4.4) follows. 
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